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I A  ABSTRACT  fC hiam  mm  0000000  KR  ft  t 


This  report  describes  the  formulation  and  use  of  a  directional  spectral 
wave  model  which  can  be  used  in  shallow  water.  The  model  simulates  the  trans¬ 
formation  of  a  directional  wave  spectrum  over  an  arbitrary  two-dimensional  ba¬ 
thymetry  allowing  for  the  effects  of  atmospheric  input,  bottom  friction,  bottom 
percolation,  nonlinear  wave/wave  interaction,  shoaling,  refraction,  and  equi¬ 
librium  range  energy  limitation.  The  user  has  the  option  to  include  or  exclude 
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these  source  and  sink  terms  in  any  combination  desired.  The  formulation  of 
each  of  these  terms  is  summarized. 

A  description  of  the  computer  code  is  provided  including  a  summary  of  the 
purpose  of  each  subroutine  and  the  format  for  data  input  and  output  of  re¬ 
sults.  A  sample  application  is  provided  showing  the  transformation  of  a  wave 
spectrum  from  deep  to  shallow  water  over  a  uniformly  sloping  bottom. 
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PREFACE 


This  report  describes  the  formulation  and  use  of  a  computer  model  to 
calculate  the  transformation  of  wave  spectra  in  shallow  water.  Various  pro¬ 
cesses  which  may  affect  wave  spectra  in  shallow  water  as  well  as  input  data 
necessary  to  initiate  calculations  are  discussed.  An  example  calculation  is 
provided  showing  input  data  and  output  results.  The  work  in  this  report  was 
authorized  by  the  Office,  Chief  of  Engineers  (OCE),  Coastal  Engineering  Func¬ 
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A  USER'S  GUIDE  TO  A  STEADY-STATE  SHALLOW-WATER 
DIRECTIONAL  SPECTRAL  WAVE  MODEL 


PART  I :  INTRODUCTION 

1.  The  need  to  predict  wave  conditions  in  shallow  water  has  existed  for 
many  years.  Recently,  shallow-water  wave  prediction  has  received  increased 
emphasis  in  response  to  the  need  for  more  accurate  predictions.  Reliable  cal¬ 
culations  of  deepwater  wave  climates  can  now  be  made  and  provide  input  for 
propagation  into  shallow  water.  Improved  shallow-water  estimates  will  improve 
the  economics  of  construction  and  navigation  along  the  coasts. 

2.  Calculated  and  measured  deepwater  and  shallow-water  wave  parameters 
are  derived  predominantly  from  wave  spectra  today.  In  the  past,  wave  parame¬ 
ters  were  estimated  considering  monochromatic  waves.  Increased  understanding 
of  wave  generation,  propagation,  and  dissipation  over  the  years  has  indicated 
that  the  spectral  approach  is  generally  more  representative  of  prototype  con¬ 
ditions.  The  availability  of  large,  fast  computers  now  makes  directional 
spectral  modeling  economically  practical.  The  physical  processes  which  affect 
waves  in  deep  and  shallow  water  are  generally  dependent  on  wave  frequency. 

Wave  energy  can  be  transferred  between  waves  of  different  frequencies,  and  it 
can  be  distributed  in  different  directions  as  well  as  at  different  frequen¬ 
cies.  These  factors  require  that  a  spectral  modeling  approach  be  adopted. 

3.  Spectral  modeling  evolved  in  deep  water  (Phillips  1957,  Pierson  and 
Moskowitz  1964,  Barnett  1968,  Ewing  1971  and  Hasselmann,  et  al.  1973)  and  is 
now  being  extended  into  shallow  water  (Hsiao  1977  and  Hughes  and  Jensen,  in 
preparation).  Shallow-water  modeling  increases  in  complexity  over  deep-water 
modeling  due  to  the  interaction  of  the  waves  with  the  bottom.  Such  processes 
as  bottom  friction,  percolation,  and  bottom  motion  may  now  be  important. 

Other  processes  present  in  deep  water  such  as  wave  interaction,  white  capping, 
breaking,  and  wind  generation  are  affected  in  shallow  water  since  wave  proper¬ 
ties  such  as  speed,  height,  and  length  change  as  a  function  of  depth.  The  re¬ 
gional  scale  of  a  shallow-water  model  also  will  be  generally  smaller  than  that 
of  a  deepwater  model  and  probably  adjacent  to  a  coastline.  This  size  differ¬ 
ence  may  result  in  more  variability  in  space  and  time  of  the  wind  field  in  a 
region  due  to  land  effects  which  adds  complexity  to  a  problem.  The  various 
processes  mentioned  above  are  referred  to  as  source  or  sink  terms  depending  on 


whether  they  add  or  subtract  energy,  respectively,  from  the  wave  field.  It  is 
common  to  refer  to  source  and  sink  terms  collectively  as  source  terms.  The 
representation  of  various  source  terms  is  still  an  area  of  research. 

k.  The  formulation  of  this  model,  which  evolved  from  Hsiao  (1977), 
will  be  discussed  below  with  a  description  of  the  present  formulation  of  each 
source  term.  Next,  a  description  of  the  computer  program  will  be  provided 
with  definitions  of  subroutines,  important  variables,  and  flow  of  the  calcula¬ 
tions.  Finally,  a  sample  calculation  will  be  presented  showing  the  structure 
of  the  input  data  and  the  calculated  results. 


PART  II:  MODEL  FORMULATION 


5.  The  model  is  based  on  the  assumption  that  the  wave  field  can  be 

represented  by  a  continuous  spectrum  of  variance  or  energy  over  a  range  of 

spectral  components  k„  ,  k  *  in  wave  number  space  or  in  frequency,  direction 

x  y 

(f,8)  space.  The  spectral  density  is  denoted  by  F  in  k  space  and  E  in 
(f,e)  space.  Each  component  can  lose  or  gain  energy,  and  intraspectral 
transfer  of  energy  due  to  weak  nonlinear  interactions  is  permitted.  The  radi¬ 
ative  transfer  equation  is  used  to  express  the  spectral  balance  in  shallow 
water.  In  wave  number  space,  this  equation  is 

DF_  3F  3Fdx  3F  dy,  3F  dky  ... 

Dt  ”  at  +  3x  dt  +  3y  dt  ^  3k  dt  +  3k  dt  '  ' 

x  y 

where  F=F(x,y  ,  k,  t)  and  S  represents  source  and  sink  terms. 
aF 

The  term  —  represents  the  time  rate  of  change  of  the  spectrum.  This  model 

3F 

assumes  that  a  steady  state  exists  so  that  =  0  .  The  next  two  terms  rep¬ 
resent  change  of  spectral  energy  in  space,  since  wave  energy  propagates  at 
the  group  velocity,  dx/dt  =  C_  and  dy/dt  =  C  where  C  is  the  group 

velocity  in  the  x  and  y  directions,  respectively.  The  last  two  terms  on 
the  left  hand  side  represent  change  of  spectral  energy  with  wave  number.  The 

terms  dk  /dt  and  dk  /dt  can  be  expressed,  respectively,  as  -aw/ax  and 
x  y 

-3w/ay  by  using  the  conservation  of  waves  equation 
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Vw 
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(2) 


where 

w  =  radial  wave  frequency  (w  =  2*f) 

V  =  horizontal  differential  operator  3/3x  +  a/3y 

6.  Using  the  steady-state  assumption  and  relations  above,  Equation  1 
becomes 


C  3F  C  3F 

g  -  +■  g  - 

Bx  ax  6y  ay 


3w  a_F_ 
ax  ak 
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3ui  aF 
ay  ak 
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(3) 


*  For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation 
(Appendix  F) . 


It  is  now  assumed  that  the  wave  frequency  does  not  change  as  a  function  of 
space  along  the  path  of  propagation,  so  Equation  3  becomes 


C  3F  C  3F  .  .... 

3x  +  3y 

The  left  hand  side  represents  refraction  and  shoaling,  and  the  right  hand  side 
represents  various  source  or  sink  terms-  This  equation  can  be  expressed  in 
frequency  direction  space  by  using,  the  following  relation: 

E(f,0)  df  =  F(ikKe)  0  (d fk| )  (5) 

=  F(k)  (k|  (dfk|  ) 

First,  express  k  in  polar  coordinates  where  |k|  is  the  radius  in  polar  co¬ 
ordinates  and  f  is  radial  frequency  to  get 

^  n  E(w,0)  dtj  d0  =  F(k,0)k  dk  d0 

Then,  integrate  over  e  from  -tt  to  t  to  get  E(u>)du»  =  2irF(k)k  dk  and  note 
that  du/dk  is  the  group  velocity  Cg  .  Thus  E  =  (2Trk/C  )F  ,  and  Equation  4 


becomes 


3x  [  k  J  8y  ay  [  k  J 


which  is  valid  for  both  deep  and  shallow  water.  This  is  the  equation  which  is 
integrated  numerically  in  this  model  to  propagate  energy  from  an  offshore 
boundary  to  the  shore. 

Propagation 

7.  If  the  source  terms  in  Equation  6  are  set  to  zero,  the  solution  rep¬ 
resents  propagation  with  refraction  and  shoaling  of  a  wave  component  along  a 
wave  ray.  In  the  absence  of  currents,  the  property 


E(f,e)  rf-  =  CONSTANT  =  E(F,0  )  ~ 


holds  where  the  zero  subscript  refers  to  initial  conditions.  This  equation  is 
solved  along  ray  paths  given  by 
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dx 

ds 


tt  =  cos  0  ,  ^  =  sin  0 


d*  a 
ds 


where  s  is  a  coordinate  along  the  ray  and 


ds  r 
dt  '  L 


contours  can  be  used,  k  sin  0  =  k  sin  0^ 

o  o 


can  be  calculated  from  this  using  k  cos  0A0  =  k  cos  0  A0 

o  o  o 


where  C  is  the  phase  velocity.  The  solution  is  obtained  row  by  row  on  the 
numerical  grid  moving  from  the  outer  boundary  toward  the  shore.  The  depth 
contours  are  assumed  to  be  locally  parallel  so  that  Snell's  Law  for  parallel 

The  change  in  wave  angle  ae 

A  limitation 

on  the  value  of  cos  0  is  imposed  so  that  singularities  do  not  occur.  The 
refraction  and  shoaling  method  in  the  model,  based  on  Snell's  Law,  is  subject 
to  the  problems  which  have  been  evident  in  monochromatic  refraction  programs. 
The  user  should  be  cautious  in  applying  the  model  to  highly  irregular  bathyme¬ 
try.  Errors  possibly  introduced  by  refraction  in  these  cases  may  not  be  evi¬ 
dent  due  to  the  limitations  imposed  and  the  directional  spectral  form  of  the 
results. 


Source  Terms 


8.  It  is  evident  from  Equation  7  that  energy  is  conserved  in  the  ab¬ 
sence  of  any  source  terms.  It  may  be  redistributed  in  (f,0)  space,  but 
nothing  is  lost  or  gained.  The  addition  of  source  terms  allows  energy  to  be 
lost,  gained,  or  redistributed  within  the  spectrum.  The  source  terms 
available  in  this  model  are  energy  input  from  winds,  energy  loss  from  bottom 
friction,  percolation,  white  capping,  and  breaking  and  redistribution  of 
energy  due  to  nonlinear  interactions.  Each  of  these  is  discussed  briefly 
below. 

Atmospheric  input 

9.  Two  mechanisms  have  been  suggested  as  responsible  for  energy  input 
to  a  wave  field  from  the  atmosphere.  Phillips  (1957)  suggested  that  turbulent 
pressure  fluctuations  in  the  atmosphere  lead  to  linear  growth  of  wave  energy. 


Thi3  is  represented  in  the  model  by  SWp  which  is  given  by 


S  =  a  t"T 
wp  g  2 nk 


(8) 
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u 
u 

=  difference  between  wind  and  wave  angles 
=  wind  speed 

The  linear  growth  coefficient  a  is  based  on  the  empirical  work  of  Snyder 

8 

(1974).  This  term  has  been  important  in  deepwater  modeling  as  a  mechanism  to 
initiate  wave  growth.  The  initial  waves  would  then  grow  through  another  me¬ 
chanism.  A  spectrum  is  input  as  an  initial  condition  in  this  model,  so  this 
term  is  of  minor  importance.  Energy  from  the  atmosphere  is  input  mainly 
through  the  following  source  term. 

10.  Miles  (1957)  suggested  a  growth  mechanism  due  to  a  coupling  between 
existing  waves  with  their  induced  pressure  fluctuations  and  the  wind.  He  hy¬ 
pothesized  that  this  would  lead  to  exponential  growth  of  waves  in  time  and  a 
shift  of  the  spectral  peak  to  lower  frequencies.  Hsiao  and  Shemdin  (1983) 
modified  this  formulation  based  on  data  from  the  Marine  Remote  Sensing  Experi¬ 
ment  (MARSEN).  This  wave  growth  mechanism  is  represented  in  the  model  as 
Swn,  which  is  given  by 

S  .0.12  —fe  cos  9  -  l\wE  (12) 

“  *w\c  / 

where  p  and  p  are,  respectively,  the  densities  of  air  and  water.  Note 
that  the  magnitude  of  this  source  term  is  dependent  on  the  ratio  of  wind  speed 
to  wave  speed.  Waves  in  a  spectrum  will  have  different  speeds  and  directions 

1  /p 

and  will  approach  a  shallow-water  wave  speed  C  =  (gh)  at  different  loca¬ 
tions  due  to  their  period  and  the  depth.  Thus,  the  contribution  of  this 
source  term  can  be  variable  across  the  spectrum  and  model  region  for  a  con¬ 
stant  wind  speed  and  direction. 


where 

q 
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Bottom  friction 

11.  In  shallow  water  there  exists  a  boundary  layer  above  the  bed 
which  in  part  is  due  to  waves.  Energy  is  dissipated  through  friction  in  this 
boundary  layer  when  water  flows  over  the  bed.  Hasselmann  and  Collins  (1968) 
proposed  an  equation  to  represent  this  frictional  loss  in  a  spectral  wave 
model.  The  frictional  source  term  is  based  on  their  formulation  and  is  repre¬ 
sented  in  the  model  as 

C  /  ?  -  W  W  \ECa 

st  *  -2  ? — 2k  <i“i>  *  kik2  <  -p=r  ^2^  (,3) 

w  cosh  kh\  |w|  K 

where  Cp  is  a  friction  coefficient  and  W  is  water  velocity  at^the  top  of 
the  boundary  layer.  The  components  of  vectors  W  and  k  are  W1  ,  W2 
and  k1  ,  k2  ,  and  <  >  is  an  ensemble  average  for  a  number  of  waves.  The 
friction  coefficient  is  determined  from  empirical  relations  (Hsiao  and 

Shemdin  1978)  and  a  specified  sediment  diameter.  Nonwave- induced  bottom  cur¬ 
rents  can  be  specified.  If  they  are  specified  as  zero,  only  wave-induced  bot¬ 
tom  currents  are  used.  The  empirically  derived  relation  between  Reynolds  num¬ 
ber,  relative  roughness,  and  friction  coefficient  is  used  to  determine  Cf  . 
Values  of  Cf  can  range  from  0.005  to  0.5  depending  on  wave  and  bottom 
characteristics.  Provision  is  included  for  the  presence  or  absence  of  ripples 
on  the  bottom.  It  is  possible  that  during  storm  conditions  the  frictional 
losses  might  suddenly  increase  with  the  onset  of  ripple  formation;  and  then  if 
the  wave  conditions  become  more  severe,  the  losses  might  decrease  as  the 
bedforms  are  eliminated. 

Percolation 

12.  Wave  motion  in  a  fluid  above  a  porous  bed  induces  an  oscillating 
flow  within  the  porous  media  which  results  in  dissipation  of  wave  energy. 
Shemdin,  et  al.  (1978)  developed  a  spectral  source  term  for  this  process  based 
on  Putnam's  (1949)  work.  This  source  term  is  represented  in  the  model  as 


S  =  -  ka 
P 


tanh  kd 
^  cosh^kh 


where 


dp  =  permeability  coefficient 
d  =  depth  of  the  permeable  layer 

It  is  assumed  that  the  permeable  layer  is  isotropic  and  of  constant  depth. 


H 
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Values  of  a  range  from  .01  to  .001  m/sec  for  coarse  to  fine  sand.  There  is 
not  a  consensus  among  shallow-water  wave  modelers  on  the  importance  of  bottom 
effects  such  as  friction,  percolation,  or  bottom  motion  in  dissipating  wave 
energy.  Studies  to  quantify  these  effects,  especially  for  irregular  waves  in 
prototype  conditions,  are  difficult  to  carry  out.  Aranuvachapun  and  Thornton 
(1984)  provide  some  insight  in  this  area  using  the  Atlantic  Remote  Sensing 
Land-Ocean  Experiment  (ARSLOE)  data  set. 

Nonlinear  wave-wave  energy  transfer 

13.  Surface  gravity  waves  can  be  described  mathematically  by  a  linear 
equation.  This  is  a  first  order  approximation  describing  some  of  the  behavior 
of  prototype  waves  but  is  not  a  complete  description.  In  a  linear  approxima¬ 
tion,  waves  travel  over  the  water  surface  independent  of  each  other  and  un¬ 
changed  in  space  and  time  assuming  that  the  depth  or  some  source  term  does  not 
affect  the  waves.  If  one  extends  the  mathematical  description  to  second  or¬ 
der,  a  transfer  of  energy  to  the  second  harmonic  occurs;  and  the  surface  pro¬ 
file  of  the  wave  changes  from  sinusoidal  to  trochoidal  form.  An  extension  to 
third  order  causes  a  change  in  phase  velocity.  It  is  apparent  from  these 
higher  order  nonlinear  descriptions  that  waves  will  interact  with  each  other, 
changing  the  waves  present  and  possibly  generating  new  waves.  These  effects 
are  termed  nonlinear  wave-wave  interactions.  They  can  be  strong  or  weak  and 
occur  in  deep  or  shallow  water. 

14.  Possibly  the  most  striking  nonlinear  interaction  is  wave  breaking 
which  can  transfer  energy  into  higher  and  possibly  lower  frequencies  than 
those  of  the  breaking  wave.  Nonlinear  interactions  have  been  successfully 
used  to  explain  the  evolution  and  shape  of  the  spectra  calculated  from  the 
Joint  North  Sea  Wave  Program  (JONSWAP)  data  (Hasselmann  1961,  1963a,  1963b). 
Third  order  nonlinear  wave-wave  interactions  of  the  type  which  contributed  to 
the  shape  of  JONSWAP  spectra  can  be  calculated,  but  the  calculations  are  so 
complex  and  time  consuming  that  it  is  not  practical  to  incorporate  them  into  a 
spectral  wave  model.  For  this  reason,  the  nonlinear  source  term  is  parameter¬ 
ized  as  a  function  of  frequency  and  spectral  width  parameters  (Hasselmann  et 
al.  1980).  This  is  applicable  in  deep  water.  In  finite  depth  water,  a  scal¬ 
ing  of  deepwater  values  is  used  based  on  computations  of  Herterich  et  al. 
(1980).  The  nonlinear  source  term  is  sensitive  to  the  energy  level  and  shape 
of  the  spectrum.  Since  both  of  these  can  be  expected  to  change  in  shoaling 
water,  a  parameterization  based  on  a  JONSWAP  spectrum  may  not  be  the  most 
appropriate. 
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Shallow-water  equilibrium 
range  limitation _ 


15.  Studies  of  waves  in  deep  water  have  led  to  the  concept  of  an 
equilibrium  range  in  the  spectrum  where  energy  input  is  in  equilibrium  with 
energy  dissipation.  Phillips  (1958)  proposed  a  spectral  form  in  this  range  as 

Ep(f)  =  ag2f“5  (2*)-11 

Pierson  and  Moskowitz  (1964)  modified  this  form  to  allow  for  a  forward  face  on 
the  spectrum  for  frequencies  lower  than  the  peak  frequency  fm  .  Their  formu¬ 
lation  is 

-5/Hr/t  r4 

EPM(f)  *  Ep<r)  e 


Data  from  JONSWAP  were  used  by  Hasselman,  et  al.  (1973)  to  propose  a  spectral 
shape  under  developing  wave  conditions  which  reaches  an  equilibrium  shape 
under  sufficient  wind  speed  and  fetch  conditions.  This  is  termed  the  JONSWAP 


spectrum  and  is  expressed  as 


Ej(f )  =  EpM(f)  Y 


-(f/f  -1)2/2o2 
e  m 


where  parameters  a  ,  y  ,  o  ,  and  fm  are  functions  of  wind  speed  and  fetch 
and 

fo  for  f  i  f 
a  m 

a  ■ 

a.  for  f  <  f 
b  m 

16.  Kitaigorodskii,  et  al.  (1975)  proposed  that  a  similar  equilibrium 
range  might  exist  in  finite  depth  water  and  proposed  a  frequency  and  depth 
dependent  factor  <j>  (w,  h)  which  would  allow  transformation  of  a  deepwater 
equilibrium  spectrum  to  one  in  finite  depth  water. 

17.  Bouws,  et  al.  (1985)  used  data  from  three  field  experiments  (Texel, 
MARSEN,  and  ARSLOE  to  propose  the  following  form  for  a  wind  sea  spectral 
shape  in  finite  depth: 


This  form,  known  as  the  TMA  spectrum,  was  fitted  to  over  2,800  wind  sea  spec¬ 
tra  from  the  above  experiments  by  varying  a  and  y  .  They  found  that  o 
and  y  could  be  expressed  empirically  as  functions  of  wind  speed  and  peak 
wave  number.  All  of  these  representations  provide  spectra  which  have  an  upper 
limit  energy  content  and  a  certain  shape  based  on  empirical  data.  The  equi¬ 
librium  range  limitation  in  the  present  model  requires  that 

E(f)  <  Ep  { f , 0 )  *  (w,h) 

This  condition  is  imposed  prior  to  the  source  terms  acting  and  represents  an 
upper  limit  of  energy  as  a  function  of  frequency  and  depth.  It  can  be  con¬ 
sidered  as  a  dissipative  source  term  where  energy  loss  is  through  turbulent 
and  viscous  processes  associated  with  white  capping  and  large  scale  breaking. 

Summary  Remarks  on  Model  Formulation 

18.  This  model  is  based  on  physical  principles,  empirical  data,  and 
assumptions  which  simplify  the  problem  enough  to  make  it  solvable.  Those 
processes  which  are  considered  most  important  in  the  propagation  and  transfor¬ 
mation  of  waves  in  shallow  water  are  included  in  the  model  formulation.  It 
should  be  recognized  that  the  representation  of  these  processes  is  the  best 
available  we  have  to  date  but  should  not  be  considered  exact  or  permanent. 

The  important  fact  that  waves  are  limited  in  energy  and  frequency  distribution 
in  shallow  water  is  represented  in  the  model.  Based  on  present  knowledge, 
wave  spectra  in  shallow  water  can  be  expected  to  remain  at  or  below  this  limi¬ 
tation.  The  relative  contributions  of  various  source  terms  used  in  shallow- 
water  wave  modeling  is  not  fully  understood.  Use  of  shallow-water  models  re¬ 
quires  care  because  the  balance  between  various  source  terms  can  vary  with  the 
situation  modeled.  Shallow-water  wave  modeling  is  not  yet  at  the  point  where 
a  model  alone  can  be  applied  at  a  site  and  be  expected  to  always  give  accurate 
wave  information  in  a  region.  It  is  recommended,  in  studies  where  these 
models  are  used,  that  wave  data  be  measured  at  various  sites  in  the  region. 
Such  data  can  be  used  to  calibrate  and  verify  the  model.  Then  the  model  can 
be  used  as  an  interpolation/extrapolation  device  and  provide  more  confidence 
that  the  results  are  correct.  More  field  measurements  will  add  to  the  data 
base  of  shallow-water  wave  measurements  and  allow  improvement  in  the  empiri¬ 
cally  based  parts  of  the  model  as  well. 


PART  III:  COMPUTER  PROGRAM 


19.  Following  is  a  description  of  the  computer  code  which  represents 
the  model  formulation  discussed  above.  The  description  applies  to  a  version 
of  the  code  which  is  operating  on  the  CDC  CYBER  200  Series  Computer  System 
compiled  with  the  FORTRAN  compiler  and  executed  with  job  control  language  (Ap 
pendix  A) .  This  code  and  computer  system  are  available  to  any  US  Army  Corps 
of  Engineers  user. 


Program  Structure 

20.  The  program  is  separated  into  two  parts.  The  first  part  calcu¬ 
lates  the  propagation,  refraction,  and  shoaling  of  all  specified  frequency  and 
direction  components  for  the  specified  bathymetry.  The  results  of  these  cal¬ 
culations  are  saved  and  used  in  the  second  part  of  the  program.  As  long  as 
the  bathymetry  and  frequency /direct ion  band  widths  remain  the  same,  the  re¬ 
sults  of  the  first  part  can  be  used  to  transform  various  spectra  over  the 
bathymetry  using  various  combinations  of  source  terms.  The  calculations  in 
the  first  part  of  the  program  are  the  most  costly  so  that  making  them  only 
once  increases  efficiency  when  making  multiple  spectral  calculations. 

21.  The  second  part  of  the  program  transforms  the  refracted  and  shoaled 

spectrum  using  one  or  more  of  the  following  source  terms: 

a.  Bottom  friction. 

b.  Percolation. 

c.  Atmospheric  input. 

(1)  Miles'  mechanism. 

(2)  Phillips'  mechanism. 

d.  Nonlinear  interactions. 

e.  Depth  limitation. 

Structuring  the  program  this  way  allows  one  to  include  those  source  terms 
which  are  most  applicable  for  a  given  problem. 

22.  The  first  part  of  the  program  consists  of  a  main  program  and  one 
subroutine  to  calculate  wave  number.  The  second  part  of  the  program  consists 
of  a  main  program  and  seventeen  subroutines.  Both  parts  have  parameter  state¬ 
ments  in  the  main  program  and  subroutines  which  allow  dimensions  of  arrays  to 
be  set.  The  variables  in  the  parameter  statement  are  defined  below. 


Definition 


Parameter  Marne 
NTX 
NTV 
NTF 
NTA 
MOUT 


Number  of  grid  cells  in  the  x  direction 
Number  of  grid  cells  in  the  y  direction 
Number  of  frequency  bands 
Number  of  angle  bands 

Number  of  grid  cells  where  printout  of  spectra  are  desired 
The  purpose  of  each  of  the  subroutines  is  described  briefly  below. 


Subroutine  Descriptions 


TWV1 

23.  TWV1  reads  in  the  parameters  to  specify  which  source  terms  are  to 
be  used.  It  calls  subroutine  TBND  and  calculates  some  constants.  It  is 
called  once  from  the  main  program. 

TBND 

24.  TBND  reads  in  the  type  of  boundary  condition  used  as  well  as  the 
directional  spectra  on  the  input  boundary.  It  is  called  once  from  TWV1 . 

TLIM 

25. TLIM  calculates  the  upper  bound  of  a  wind  sea  spectrum  as  a  function  of 
depth  and  frequency.  It  is  called  at  each  grid  point  from  the  main  program. 

TSPF 

26.  TSPF  locates  the  frequency  in  the  spectrum  at  which  the  maximum 
energy  occurs  and  calculates  the  significant  wave  height  based  on  the  total 
energy  in  this  spectrum.  It  is  called  at  each  grid  point  and  in  turn  calls 
TCFT. 

TSPA 

27.  TSPA  calculates  the  angle  at  which  the  maximum  energy  in  the  spec¬ 
trum  occurs.  It  is  called  at  each  grid  point  and  in  turn  calls  TCFT. 

TCFT 

28.  TCFT  calculates  a  quadratic  fit  to  part  of  the  energy  spectrum  to 
determine  where  the  peak  frequency  and  direction  occur.  It  is  called  from 
TSPF  and  TSPA. 

WAVENO 

29.  WAVENO  calculates  the  wave  number,  given  the  frequency  and  depth, 
using  the  dispersion  relation  for  linear  waves.  It  is  called  in  the  main  pro¬ 
gram  and  TFRC,  TNL,  and  TLIM. 


TDIS 


30.  TDIS  is  the  first  subroutine  in  the  calculation  of  bottom  friction. 
It  calculates  the  ratios  of  bottom  velocity  components  used  in  the  formulation 
of  bottom  friction.  It  is  called  from  the  main  program  at  every  grid  point. 
TFRC 

31.  TFRC  is  the  second  subroutine  in  the  bottom  friction  calculation. 

It  calculates  the  bottom  friction  coefficient  based  on  empirical  data  and  a 
knowledge  of  the  sand  diameter.  The  bottom  friction  coefficient  may  be  spec¬ 
ified  directly  also.  It  is  called  from  the  main  program  at  each  grid  point. 
TNL 

3 2.  TNL  calculates  the  nonlinear  interactions  using  a  parameterization 
based  on  exact  calculations  of  the  nonlinear  interactions  for  a  JONSWAP  spec¬ 
trum.  A  similarity  relationship  allows  scaling  to  finite  depth  water  for 
kph  £  0.7  .  It  is  called  from  the  main  program  at  every  grid  point. 

TPER 

33.  TPER  calculates  the  percolation  source  term  with  a  specified  value 
of  the  permeability  coefficient  and  depth  of  the  sand  layer.  It  is  called 
from  the  main  program  at  every  grid  point,  frequency,  and  angle. 

TFRA 

34.  TFRA  is  the  third  part  of  the  bottom  friction  coefficient.  It  cal¬ 
culates  a  product  of  terms  in  the  bottom  friction  calculation.  It  is  called 
from  the  main  program  at  every  grid  point  and  every  frequency. 

TFR 

35.  TFR  is  the  final  part  of  the  bottom  friction  coefficient  in  which 
various  terms  are  combined.  It  is  called  from  the  main  program  at  every  grid 
point,  frequency,  and  angle. 

TWN 1 

36.  TWN1  is  the  calculation  of  the  Miles  atmospheric  energy  transfer 
mechanism.  It  is  called  from  the  main  program  at  every  grid  point,  frequency, 
and  angle. 

TWN  2 

37.  TWN2  is  the  calculation  of  the  Phillips  atmospheric  energy  transfer 
mechanism.  It  is  called  from  the  main  program  at  every  grid  point,  frequency, 
and  angle. 

TWN2A 

38.  TWN2A  calculates  parameters  used  in  the  formulation  of  the  Phillips 


atmospheric  energy  transfer  mechanism.  It  is  called  from  the  main  program  at 
every  grid  point  and  frequency. 

TWV3 

39.  TWV3  prints  out  arrays  of  significant  height,  peak  frequency,  and 
peak  direction  after  all  calculations  are  finished.  It  is  called  once  from 
the  main  program. 


Input  Data 

40.  There  are  two  sets  of  input  data,  one  each  for  the  refraction  and 
transformation  parts  of  the  program.  The  input  data  for  the  refraction  part 
of  the  program  consists  of  the  following: 

a.  Up  to  four  80-column  lines  of  descriptive  text. 

b.  Namelist  SIZE. 

DX  -  Grid  cell  size  in  x  direction  in  meters 
DY  -  Grid  cell  size  in  y  direction  in  meters 

DELSI  -  Incremental  distance  along  ray  path  in  meters  (typi¬ 
cally  one-fourth  the  grid  cell  size) 

c.  Namelist  ANGLE. 

THETA  -  Directions  in  degrees  of  the  center  of  angular  band3  of 
width  A6=e2-e1  (Bands  must  be  symmetrically  placed  about 
0  deg  which  is  normal  to  shore.) 

d.  Namelist  FREQU. 

FREQ  -  Center  of  frequency  bands  in  hertz 
DFR  -  Bandwidth  in  hertz  of  each  band 

e.  Namelist  DEPTH. 

HT  -  Depth  in  meters  at  each  grid  point  (A  positive  value 
indicates  depth  below  mean  sea  level.  Depths  are  read  in 
columns  left  to  right  from  offshore  to  shore.) 

HTCNV  -  A  conversion  factor  to  convert  depths  to  meters  if 
depths  read  in  are  not  in  meters. 

A  typical  grid  orientation  is  shown  in  Figure  1.  An  example  of  an  input  file 
for  this  part  of  the  program  is  shown  in  Appendix  B. 

41.  The  input  data  for  the  transformation  part  of  the  program  consist 
of  specification  of  the  types  of  source  terms  and  boundary  conditions  desired; 
the  values  of  spectral  energy  as  a  function  of  frequency  and  direction  speci¬ 
fied  on  the  input  boundary;  and  three  namelists  supplying  coefficients  for  the 
frictional,  percolation,  and  atmospheric  input  subroutines  and  indices  of  grid 


V 


SHORELINE 

Figure  1 .  Orientation  of  model  grid  with 
convention  for  wave  direction 


points  at  which  spectral  output  is  desired. 

42.  A  source  term  can  be  turned  on  or  off  by  specifying  T  or  F  , 
respectively,  for  true  or  false.  The  format  is  (LI,  5  (IX,  L2)),  and  in  se¬ 
quential  order  the  source  terms  are  friction,  percolation,  Mile's  atmospheric 
input,  Phillips'  atmospheric  input,  nonlinear  interactions,  and  depth  limita¬ 
tion  . 


43.  The  types  of  boundary  conditions  are  specified  by  the  values  of 
three  parameters  IDS,  IDL,  and  IDU  which  represent,  respectively,  the  side, 
lower,  and  upper  boundaries  as  shown  in  Figure  1.  The  following  values  of 
IDS,  or  IDL,  or  IDU  will  result  in  the  following  action: 

a.  The  energy  value  of  all  spectral  components  (frequency  and  di¬ 
rection)  at  each  grid  point  on  the  specified  boundary  will  be 
set  equal  to  a  constant  specified  with  the  next  read  statement 
using  a  G 1 3 - 5  format. 

b.  The  spectrum  is  equal  at  each  point  along  the  boundary  and 
specified  with  the  next  read  statement  using  a  6G13.5  format. 
Values  are  read  from  lowest  to  highest  frequency  at  a  given  an¬ 
gle  starting  with  the  smallest  angle,  where  the  smallest  angle 
is  the  largest  negative  value. 

c.  The  spectrum  is  read  in  at  each  point  on  the  particular  bound¬ 
ary  using  6G13.5  format. 

d.  The  spectral  values  on  the  upper  and  lower  boundaries  are  set 
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equal  to  the  values  at  the  adjacent  grid  points.  No  values  are 
read  in. 


44.  If  printout  of  one-  and  two-dimensional  spectra  are  desired,  enter 
the  x-  and  y-coordinate  indices  of  the  point  in  Format  (IX,  213).  The  final 
data  input  is  through  three  namelists  called  FRICT,  PERCO,  AND  GENER. 

a.  Namelist  FRICT. 

DSC  -  Sand  grain  diameter  used  in  the  calculation  of  the 
bottom  friction  coefficient  in  meters 

UM1  -  The  x  component  of  steady  nonwave  induced  bottom  cur¬ 
rent  in  meters  per  second 

UM2  -  The  y  component  of  steady  nonwave  induced  bottom  cur¬ 
rent  in  meters  per  second 

CF1  -  A  specified  value  of  bottom  friction  coefficient  (If  a 
value  of  zero  is  specified,  the  program  will  calculate  an  ap¬ 
propriate  value.) 

b.  Namelist  PERCO. 

PKC  -  The  permeability  coefficient 

DC  -  The  depth  of  the  permeable  layer  in  meters 

c.  Namelist  GENER. 

THWC  -  The  direction  toward  which  the  wind  is  blowing  in  de¬ 
grees  with  the  same  convention  as  the  wave  approach  angle 

WC  -  The  wind  speed  in  meters  per  second 

An  example  of  an  input  file  for  this  part  of  the  program  is  shown  in  Appen¬ 
dix  C. 

Output 

45.  The  output  of  the  refraction  part  of  the  program  consists  of  line 
printer  and  disk  file  output  for  subsequent  use  by  the  transformation  part  of 
the  program.  The  following  information  is  printed  out:  the  date  and  time  the 
model  was  run;  descriptive  text  describing  the  run;  content  of  the  three  name- 
lists — SIZE,  ANGLE,  FREQU;  the  depths  at  each  grid  point;  and  a  conversion 
factor  to  meters  if  the  depths  are  not  input  in  meters.  An  example  of  printed 
output  is  shown  in  Appendix  D.  The  user  does  not  normally  interact  with  the 
disk  file  output;  consequently,  it  is  not  described  in  detail. 

46.  The  output  of  the  transformation  part  of  the  program  consists  of 
the  following  printout:  the  data  and  time  the  model  was  run,  descriptive 
text,  the  water  depths  at  each  grid  point,  the  true  or  false  status  of  each  of 


the  source  terras,  index  values  for  the  types  of  boundary  conditions,  the  spec¬ 
tral  values  on  the  boundaries,  the  significant  wave  height  at  each  grid  point, 
and  the  peak  frequency  and  peak  direction  at  each  grid  point.  An  example  of 
this  output  is  shown  in  Appendix  E.  Additional  output  of  one-  and  two- 
dimensional  spectra  at  specified  grid  points  is  an  option.  This  is  initi¬ 
ated  in  the  parameter  statement  of  the  program  by  specifying  a  value  for 
MOUT,  the  number  of  grid  points  where  output  is  desired.  The  locations 
of  the  grid  points  are  specified  after  the  spectral  energy  input.  This  is 
shown  in  Appendix  C.  An  example  of  the  output  is  shown  in  Appendix  E. 


PART  IV:  DESCRIPTION  OF  SAMPLE  CALCULATIONS 


47.  The  sample  results  are  shown  to  provide  the  user  some  feeling  for 
the  effects  of  the  source  terms  and  various  values  of  coefficients.  A  simple 
bathymetry  was  chosen  which  consists  of  a  plane  beach  with  maximum  offshore 
depth  of  50  m  and  minimum  depth  at  the  shore  of  2.5  m.  The  bathymetry  was 
approximated  with  a  20  by  20  grid.  The  grid  cells  were  of  equal  dimension 
(500  m)  in  the  x  and  y  directions.  Twelve  angular  bands  of  width  15  deg  and 
symmetric  around  0  deg  were  specified.  Twenty  frequency  bands  of  width  0.01 
Hz  from  0.06  to  0.25  Hz  were  specified.  A  J0NSWAP  spectrum  with  the  following 
parameters  was  input  on  the  offshore  boundary  with  the  mean  wave  direction 
normal  to  shore:  a  =  0.01  ,  y  =  3.3  ,  fm  =  0.1  ,  aa  =  0*07,  and  =  0.09  . 

A  J0NSWAP  spreading  function  was  used  to  distribute  energy  over  the  angular 
bands.  This  spectrum  was  applied  equally  at  each  grid  point  on  the  offshore 
boundary;  hence  IDS  =  2  .  The  upper  and  lower  boundary  conditions  were 
supplied  from  adjacent  grid  points  by  setting  IDL  =  IDU  =  4  . 

48.  Results  are  summarized  by  plots  of  significant  wave  height  at  each 

grid  point  on  a  line  perpendicular  to  shore  and  midway  in  the  grid  in  the 
alongshore  direction.  Figures  2-9  show  the  effect  of  each  source  term  acting 
alone  as  well  as  no  source  terms  acting.  Refraction  and  shoaling  are  included 

in  all  cases.  A  quantitative  summary  is  provided  in  Table  1.  Two  things  are 

obvious  from  these  figures.  First,  there  is  little  difference  between  the  ef¬ 
fect  of  each  source  term  with  exception  of  the  depth  limitation  term.  The  ef¬ 
fect  of  each  source  term  is  dependent  on  either  a  coefficient,  wind  velocity, 
or  depth.  For  the  coefficients  and  wind  value  chosen  in  this  case,  there  is 
little  difference  over  refraction  and  shoaling  alone.  Choice  of  different  co¬ 
efficient  values  can  result  in  more  significant  effects.  The  second  charac¬ 
teristic  to  note  is  the  unrealistic  behavior  of  the  wave  heights  in  depths 
where  they  would  be  expected  to  break  if  a  breaking  mechanism  is  not  included. 

49.  In  an  application  of  the  model,  it  is  important  to  have  prototype 

measurements  so  the  model  can  be  calibrated  by  adjusting  coefficient  values. 

In  the  absence  of  any  prototype  measurements,  it  is  probably  best  to  operate 
the  model  in  the  most  conservative  mode  and  employ  only  the  depth  limitation 
with  refraction  and  shoaling.  If  one  has  knowledge  of  bottom  characteristics 
or  atmospheric  input,  the  appropriate  source  terms  can  be  applied.  Figure  10 
provides  an  indication  of  the  effect  of  an  upper  limit  frictional  coefficient 
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Figure  8.  Wave  height  normal  to  shore  with 
only  depth  limitation  acting 


Figure  9.  Wave  height  normal  to  shore  with  all 
source  and  sink  terms  acting 


Significant  Wave  Height*  as  a  Function  of  Location  and  Source  Terms 
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A  BOTTOM  FRICTION  ONLY  CF  -  0.5 
•  OEPTH  LIMITATION  ONLY 


Figure  10.  Wave  height  normal  to  shore  for 
different  values  of  bottom  friction  and  for 
depth  limitation  only 


on  wave  heights  for  this  sample  case.  An  order  of  magnitude  increase  in  the 
percolation  coefficient  for  this  case  has  negligible  effect.  However,  for  a 
different  case,  such  as  a  large  extent  of  flat  bottom  with  high  permeability, 
it  could  be  important. 

50.  The  state  of  the  art  of  modeling  shallow-water  waves  has  markedly 
improved  as  exemplified  by  this  model.  Calculation  of  refraction  and  shoaling 
of  monochromatic  waves  following  rays  has  evolved  to  more  sophisticated  numer¬ 
ical  routines  for  refraction  and  shoaling  of  spectra.  A  variety  of  source 
terms  representing  possible  sources  and  sinks  of  energy  have  been  added  to 
models  to  account  for  wave  generation  and  dissipation  in  shallow  water.  Re¬ 
search  continues  on  the  most  appropriate  formulation  for  these  source  terms 
and  their  balance  in  a  spectral  model.  The  users  of  models  such  as  this 
should  exercise  care  to  obtain  adequate  calibration  and  verification  of  the 
model,  understand  what  source  terms  dominate,  and  use  common  sense  in  inter¬ 
pretation  of  the  results. 
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APPENDIX  A:  SAMPLE  JOB  CONTROL  LANGUAGE 


FRONT  END  JCL 


CYBER  200  JCL 


/JOB 

/NOSEQ 

PCJOB ,T800 ,P4. 

/USER 

/CHARGE 

GET .TMANEW2 , PCDATN1 , PCNEWR2 
SK I  PR , PCNEWR2 . 

COPYBR ,TMANEW2 , PCNEWR2 , 1 . 
COPYBR , PCDATN1 , PCNEWR2 , 1 . 
SUBM I T , PCNEWR2 , T . 

END  OF  FILE 


PNWJOB. 

USER(U=  , PA=  )ADY 

RESOURCE( JCAT =P4 ,TL=800 ) 
CHARGE(  ,  ) 

PATTACH , LINKPOOL . 

PURGE, NEW55. 

COPYL, INPUT, T1 ,PART=G,NUM=1 . 
COPYL , INPUT , D1 , PART=G ,NUM=1 . 
FORTRAN , I =T1 , B=BT/#300 , 0=B0U . 
PATTACH, FTNUTIL. 

LOAD , BT , CN*TMAGO/#400 , GRLP=*A 
REQUEST, NEW55/8000 ,RT=W. 
TMAG0(TAPE5=D1 ,TAPE55=NEW55) 
DEFINE .NEW55/8000 ,RT=W. 
RETURN, T1 ,D1 . 

END  OF  FILE 


Job  Control  Language  to  Run  the  Transformation  Program  Onlj 


FRONT  END  JCL 

/JOB 

/NOSEQ 

PETEJB,T700,P4. 

/USER* 

/CHARGE 

GET,TAPE6. 

RENAME , PCDATA-TAPE6 . 

GET ,TMANEW2 . 

SK I  PR , TMANEW2 ,  1 ,  C . 

SKIPR, PCDATA, 1 ,C. 

GET , PNW2052 . 

SKI  PR , PNW2052 . 

COPYBR , TMANEW2 , PNW2052 , 1 
COPYBR , PCDATA , PNW20  32,1. 
SUBMIT ,PNW2052,T . 

END  OF  FILE 


CYBER  200  JCL 

TMJ . 

USER(U»  ,PA=  ) ADY 

RESOURCE( JCAT-P4 ,TL=600 ) 
CHARGEC  ,  ) 

PATTACH , LINKPOOL . 

ATTACH, TAPE35. 

COPYL , INPUT ,T , PART =G ,NUM=1 . 
COPYL, INPUT, D,PART=G,NUM=1 . 
FORTRAN , I *T , B*BT/#300 , 0=B0VLK , 
LOAD , BT , CN*TMAGO/#400 . 

TMAGO( TAPES =D) 

RETURN, T,D. 

END  OF  FILE 


APPENDIX  B:  SAMPLE  INPUT  FOR 
REFRACTION  RUN 


-22. 5, -7. 5.  •>.5,22.5,37.5,52.5,67.5,82.5 


0.12,0.13,0.14,0.15, 

.22,0.23,0.24,0.25, 


SENSITIVITY  OF  TMANEW 
JONSWAP  SPECTRUM  INPUT 
1/200  PLANE  BEACH 
DECEMBER  1984 

&SIZE  DX=500 . 0 , DY=500 .0 ,DELSI=125.0  SEND 
& ANGLE 

THETA=-82 . 5 ,  -67 . 5 , -52 . 5 , -37 . 5 , -22 . 5 , -7 . 5 .  5 , 22 . 5 , 37 . 5 , 52 . 5 , 67 . 5 ,  i 

SEND 

&FREQU  FREQ=0 .06,0 .07,0 .08,0.09,0.10,0.11,0.12,0.13,0.14,0.15, 
0.16,0.17,0.18,0.19.0.20,0.21,0.22,0.23,0.24,0.25, 
DFR=20*0 .01  SEND 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 

50.000  47.500  45.000  42.500  40.000  37.500  35.000  32.500  30.000  27. 

25.000  22.500  20.000  17.500  15.000  12.500  10.000  7.500  5.000  2. 
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.20985 

.82864 

.39877 

.17770 

.28658 

.36761 

.14997 

3.3716 

.77000 

.24308 

1 . 3849 
.65090 
.19892 

9.6882 

1.2631 

.30768 

4.5293 
1 .0003 
.24601 

20.656 

1.8014 

.36434 

10.596 

1.3614 

.28651 

33.707 

2.2663 

.40643 

18.360 

1.6618 

.31619 

42.874 

2.5369 

.42886 

24.053 

1.8328 

.33188 

42.874 

2.5369 

.42886 

24.053 
1 .8328 
.33188 

33.707 

2.2663 

.40643 

1 8 . 360 
1.6618 
.31619 

20.656 

1 .3014 
.36434 

10.596 

1.3614 

.28651 

9.6882 

1.2631 

.30763 

4.5293 
1 .0003 
.24601 

3.3716 
. 77000 
.24308 

1 . 3849 
.65090 
.19892 

.32864 

.39877 

.17770 

.23658 

.36761 

.14997 

20  10 

&FRICT  OSC= .0002,UM1=0 .0 ,UM2=0 .0  ,CF1=0 
&PERCO  PKC=.010 ,DC=5.0  6END 
&GENER  WC=15.0 ,THWC=0 .0  SEND 
END  OP  FILE 
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10941  REFRACTION  MODEL  MM  AT  11/1V84  11.08.37 


MITlVmr  OF  DMNBi 
JMMRF  SFECTWH  UfUT 
1/2M  PUOC  BEACH 
NOUEMOt  1M4 


I  20  ST  20 60 ID 
NUKR  OF  FREQUENCIES  *  20 
NMER  OF  OIRECTIMS  *  12 


(SIZE 

OK  ■  .5E+03, 

(nr  *  .3E«3, 

DELSI  *  .123E+03, 

(END 


THETA  *  -.823E+02,  -.S73E+02,  -.S25E+02,  -.373E+02,  -.225E+02,  -.7SE+01,  .75E+01,  .225E+02,  .373E+02,  .525E+02,  .673E+ 

02, 

.825E+02, 

(END 

1 

(FREQU 

FREQ  «  .SE-01,  .7E-01,  .8E-01,  .9E-01,  .1E+00,  .11E+00,  .12E+00,  .13£«Q,  .14E+00,  .15E+00,  .1SE+00,  .17E+00,  .1SEF 


.19E+00,  . 2E+00 , 
DFR  *  .  IE-01, 
-01, 

.IE-01,  .IE-.., 

(EM) 

30.000  47.500  45.000 
25.000  22.500  20.000 
50.000  47.500  45.000 
25.000  22.500  20.000 
30.000  47.500  43.000 
23.001  22.500  20.000 
30.000  47.500  45.000 
23.000  22.500  20.000 
30.000  47.500  45.000 
23.000  22.500  20.000 
50.000  47.500  45.000 
25.000  22.500  20.000 
50.000  47.500  45.000 
25.000  22.500  20.000 
50.000  47.500  45.000 
25.000  22.500  20.000 
50.000  47.500  45.000 
25.000  22.500  20.000 
30.000  47.500  45.000 
25.000  22.500  20.000 
50.000  47.500  45.000 
25.000  22.500  20.000 
50.000  47.300  43.000 


.21E+00,  .22E+00,  . 23E+00 ,  .24E+00,  .2SE400, 

.IE-01,  .IE-01,  .IE-01,  . IE-01 ,  .IE-01,  .IE-01,  .IE-01,  .IE-01,  .IE-01,  .IE-01,  .IE-01,  .IE-01,  .IE 


IE-01,  .IE-01,  .IE-01,  .IE-01, 


V.V\MVl%V, 


nniiivvffveiwsTVfrvTw 


rt.uou 

22.900 

29.UOO 

17.900 

13.000 

12.300 

10.000 

7.500 

S.000 

2. SOU 

30.009 

47.300 

45.000 

42.300 

40.000 

37.300 

33.000 

32.300 

30.000 

27.300 

23.000 

22.300 

20.000 

17.500 

15.000 

12.300 

10.000 

7.500 

5.000 

2.500 

30.000 

47.300 

43.000 

42.300 

40.000 

37.300 

33.000 

32.500 

39.000 

27.500 

23.000 

22.500 

20.000 

17.500 

15.000 

12.300 

10.000 

7.300 

3.000 

2.300 

30.000 

47.300 

43.000 

42.300 

40.000 

37.300 

33.000 

32.500 

30.000 

27.300 

23.000 

22.300 

20.000 

17.300 

13.000 

12.300 

10.000 

7.500 

3.000 

2.300 

30.000 

47.300 

43.000 

42.300 

40.000 

37.300 

33.000 

32.300 

30.000 

27.300 

23.000 

22.300 

20.000 

17.300 

15.000 

12.300 

10.000 

7.500 

5.000 

2.300 

30.000 

47.300 

43.000 

42.300 

40.000 

37.300 

33.000 

32.500 

30.000 

27.300 

23.000 

22.300 

20.000 

17.300 

15.000 

12.300 

10.000 

7.300 

3.000 

2.300 

30.000 

47.300 

43.000 

42.300 

40.000 

37.500 

33.000 

32.300 

30.000 

27.300 

23.000 

22.300 

20.000 

17.500 

15.000 

12.300 

10.000 

7.300 

5.000 

2.500 

30.000 

47.300 

43.000 

42.300 

40.000 

37.300 

33.000 

32.500 

30.000 

27.300 

25.000 

22.300 

20.000 

17.300 

15.000 

12.500 

10.000 

7.500 

3.000 

2.300 

30.000 

47.300 

45.000 

42.300 

40.000 

37.500 

33.000 

32.300 

30.000 

27.300 

23.000 

22.300 

20.000 

17.300 

13.000 

12.300 

10.000 

7.300 

5.000 

2.500 

D3 


SENSITIVITY  OF  TMANEW 
JONSWAP  SPECTRUM  INPUT 
1/200  PLANE  BEACH 
NOVEMBER  1984 

WATER  OEPTH  IN  METERS 


3  10  11  12  13  14 


20 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

19 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

18 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

17 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

16 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

15 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

14 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

13 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

12 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

11 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

10 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

9 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

3 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

7 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

6 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

5 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

4 

50.0 

47.5 

45.0 

42.5 

40  .0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

3 

50.0 

47.5 

45.0 

42.5 

40  .0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

2 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

1 

50.0 

47.5 

45.0 

42.5 

40.0 

37.5 

35.0 

32.5 

30.0 

27.5 

25.0 

22.5 

20.0 

17.5 

15.0 

WATER  DEPTH  IN 

METERS 

16  17  18  19  20 

20 

12.5 

10.0 

7.5 

5.0 

2.5 

19 

12.5 

10.0 

7.5 

5.0 

2.5 

19 

12.5 

10.0 

7.5 

5.0 

2.5 

17 

12.5 

10.0 

7.5 

5.0 

2.5 

16 

12.5 

10.0 

7.5 

5.0 

2.5 

15 

12.5 

10.0 

7.5 

5.0 

2.5 

14 

12.5 

10.0 

7.5 

5.0 

2.5 

13 

12.5 

10.0 

7.5 

5.0 

2.5 

12 

12.5 

10.0 

7.5 

5.0 

2.5 

11 

12.5 

10.0 

7.5 

5.0 

2.5 

10 

12.5 

10.0 

7.5 

5.0 

2.5 

9 

12.5 

10.0 

7.5 

5.0 

2.5 

8 

12.5 

10.0 

7.5 

5.0 

2.5 

7 

12.5 

10.0 

7.5 

5.0 

2.5 

6 

12.5 

10.0 

7.5 

5.0 

2.5 

5 

12.5 

10.0 

7.5 

5.0 

2.5 

4 

12.5 

10.0 

7.5 

5.0 

2.5 

3 

12.5 

10.0 

7.5 

5.0 

2.5 

2 

12.5 

10.0 

7.5 

5.0 

2.5 

1 

12.5 

10  .0 

7.5 

5.0 

2.5 

llV 


•J; 

V 

§ 


k 

a 


CS,  L 


IDS-  2 


OBNDS 


1 

.94223E-02 

.30739 

.89999 

1 .1047 

.32864 

.  28638 

.28494 

.34778 

.39800 

.41232 

.39877 

.36761 

.32808 

.12639 

.28640 

.10651 

.24633 

.20983 

.17770 

.14997 

2 

.11240E-01 

.42749 

1.3869 

2.7379 

3.3716 

1.3849 

1.0313 

1 .0113 

.97713 

.88630 

.77000 

.63090 

.34093 

.16323 

.44310 

.13441 

.36430 

.29733 

.24308 

.19892 

3 

.  12833E-01 

.34783 

2.4312 

3.4884 

9 . 6882 

4.3293 

2.7136 

2.2372 

1.9202 

1 .3734 

1 . 2631 

1.0003 

. 78798 
.19786 

.62012 

.16011 

.48894 

.33693 

.30768 

.24601 

4 

.  14116E-01 

. 63434 

3.3013 

8.9912 

20.636 

10.396 

3.4313 

4.0148 

3.1173 

2.3800 

1.8014 

1.3614 

1.0320 

.22713 

. 78664 
.18131 

.60383 

.46713 

.36434 

.23651 

3 

.  13012E-01 

.73437 

4.0240 

12.373 

33.707 

18.360 

8.3084 

3.8268 

4 . 2646 

3.1080 

2.2663 

1.6618 

1 .2288 
.24833 

.91747 

.19686 

.69213 

.32767 

.40643 

.31619 

6 

. 13473E-01 

.77709 

4.4330 

14.474 

42.374 

24.033 

10.608 

6 . 9968 

4.9745 

3.3434 

2.5369 

1.3323 

1 .3388 
.23949 

.98930 

.20487 

.74017 

.36022 

.42886 

.33188 

7 

. 13473E-01 

.77709 

4 . 4330 

14.474 

42.374 

24.033 

10.608 

6 . 9968 

4.9743 

3.5434 

2.5369 

1.8328 

1 . 3338 
.23949 

. 98930 
.20487 

.74017 

.36022 

.42386 

.33138 

8 

. 13012E-01 

.73437 

4.0240 

12.373 

33.707 

13.360 

8.3084 

3 . 8268 

4 . 2646 

3.1080 

2.2663 

1.6618 

1.2288 

.24833 

.91747 

.19686 

.69215 

.32767 

.40643 

.31619 

9 

. 14116E-01 

.63434 

3.3013 

3.9912 

20.656 

10.596 

3.4313 

4.0148 

3.1173 

2.3300 

1.3014 

1.3614 

1 .0320 
.22713 

. 78664 
.18131 

.60383 

. 4671 3 

.36434 

.28651 

10 

. 12833E-01 

.34783 

2.4312 

5.4384 

9. 6832 

4.5293 

2.7136 

2.2372 

1.9202 

1 . 5'’54 

1.2631 

1 .0003 

.73798 

.19786 

.62012 

.16011 

.43894 

. 38693 

.30763 

. 2460 : 

11 

. 1124QE-01 

.42749 

1.5369 

2 . 7379 

3.3716 

1.3849 

1.0313 

1.0113 

.97713 

.33630 

.77000 

.65090 

.34093 

.16323 

.44510 

.13441 

. 36430 

.29755 

.24303 

.19892 

12 

. 94223E-02 

.30739 

.39999 

1.1047 

.82364 

.28633 

.28494 

.34778 

.39800 

.41252 

.39377 

.32808 

.12639 

.28640 

.10631 

.24633 

.20933 

.17770 

.14997 

K 

I  OUT 

JOUT 

1 

1 

10 

2 

3 

10 

3 

10 

10 

4 

13 

10 

3 

20 

10 

&FRICT 

DSC 

9 

•2E-03 

UM1 

■ 

.0, 

UM2 

* 

.0, 

CF1 

■ 

. 3E+00 

SEND 

END  OF  DATA  CHECK 


RESULTS  FOR  5  OBSERVATION  POINTS  FOLLOW... 


IX  IY  DEPTH  (M)  PEAK  FREQ . < HZ)  PEAK  DIR. C DEG)  SIG.  HT.(M) 
1  10  30.00  .10  .00  5.41 

1-D  2-D 

SPCTRM  SPECTRA 


.04 

.01 

.01 

.01 

.01 

.02 

.02 

.02 

.02 

.01 

.01 

.01 

.01 

1.81 

.31 

.43 

.33 

.63 

.73 

.78 

.78 

.73 

.65 

.53 

.43 

.31 

8.73 

.90 

1.59 

2.43 

3.30 

4.02 

4.44 

4.44 

4.02 

3.30 

2.43 

1 .59 

.90 

23.66 

1.10 

2.76 

3.49 

8.99 

12.37 

14.47 

14.47 

12.37 

8.99 

5.49 

2.76 

1.10 

38.19 

.83 

3.37 

9.69 

20.66 

33.71 

42.87 

42.87 

33.71 

20.66 

9.69 

3.37 

.83 

31.00 

.29 

1.38 

4.33 

10.60 

18.36 

24.05 

24.05 

18.36 

10.60 

4.53 

1.38 

.29 

14.96 

.28 

1.03 

2.71 

3.43 

8.51 

10.61 

10.61 

8.51 

3.43 

2.71 

1.03 

.28 

10.71 

.33 

1.01 

2.26 

4 .  Okl 

3.83 

7.00 

7.00 

5.83 

4.01 

2.26 

1.01 

.35 

8.20 

.40 

.98 

1.92 

3.12 

4.26 

4.97 

4.97 

4.26 

3.12 

1.92 

.98 

.40 

6.23 

.41 

.89 

1.58 

2.38 

3.11 

3.34 

3.34 

3.11 

2.38 

1 .58 

.89 

.41 

4.73 

.40 

.77 

1.26 

1.80 

2.27 

2.54 

2.54 

2.27 

1.80 

1.26 

.77 

.40 

3.60 

.37 

.63 

1.00 

1.36 

1.66 

1.83 

1.83 

1.66 

1.36 

1.00 

.65 

.37 

2.75 

.33 

.54 

.79 

1.03 

1.23 

1.34 

1.34 

1.23 

1.03 

.79 

.54 

.33 

2.12 

.29 

.43 

.62 

.79 

.92 

.99 

.99 

.92 

.79 

.62 

.45 

.29 

1.64 

.23 

.36 

.49 

.60 

.69 

.74 

.74 

.69 

.60 

.49 

.36 

.25 

1 .28 

.21 

.30 

.39 

.47 

.53 

.56 

.56 

.53 

.47 

.39 

.30 

.21 

1.01 

.18 

.24 

.31 

.36 

.41 

.43 

.43 

.41 

.36 

.31 

.24 

.18 

.80 

.13 

.20 

.25 

.29 

.32 

.33 

.33 

.32 

.29 

.23 

.20 

.15 

.64 

.13 

.16 

.20 

.23 

.25 

.26 

.26 

.25 

.23 

.20 

.16 

.13 

.32 

.11 

.13 

.16 

.18 

.20 

.20 

.20 

.20 

.18 

.16 

.13 

.11 

SIGNIFICANT  WAVE  HEIGHTS  IN  METERS 


1  2  3  4  3  6  7  8  9  10  11  12  13  14  15 


20 

5.41 

5.29 

5.20 

5.09 

4.96 

4.82 

4.66 

4.46 

4.23 

3.98 

3.69 

3.36 

3.00 

2.62 

2.21 

19 

3.41 

3.29 

5.20 

5.09 

4.96 

4.82 

4.66 

4.46 

4.23 

3.98 

3.69 

3.36 

3.00 

2.62 

2.21 

18 

5.41 

5.29 

5.20 

5.09 

4.96 

4.32 

4.66 

4.46 

4.23 

3.98 

3.69 

3.36 

3.00 

2.62 

2.21 

17 

5.41 

5.31 

5.20 

5.09 

4.97 

4.82 

4.66 

4.46 

4.23 

3.98 

3.69 

3.36 

3.00 

2.62 

2.21 

16 

5.41 

5.31 

5.20 

3.08 

4.96 

4.82 

4.65 

4.45 

4.23 

3.97 

3.63 

3.36 

3.00 

2.61 

2.20 

15 

5.41 

5.31 

5.20 

5.09 

4.96 

4.82 

4.65 

4.43 

4.23 

3.97 

3.68 

3.36 

3.00 

2.61 

2.21 

14 

5.41 

5.31 

5.20 

5.08 

4.96 

4.32 

4.65 

4.43 

4.23 

3.97 

3.68 

3.36 

3.00 

2.61 

2.20 

13 

5.41 

5.31 

5.20 

5.09 

4.96 

4.81 

4.64 

4.43 

4.23 

3.97 

3.68 

3.36 

3.00 

2.61 

2.20 

12 

5.41 

5.32 

5.20 

5.09 

4.96 

4.81 

4.64 

4.45 

4.23 

3.97 

3.68 

3.36 

3.00 

2.61 

2.20 

11 

5.41 

5.32 

5.20 

5.09 

4.96 

4.81 

4.64 

4.45 

4.23 

3.97 

3.68 

3.36 

3.00 

2.81 

2.20 

10 

5.41 

5.32 

5.20 

5.09 

4.96 

4.81 

4.64 

4.45 

4.23 

3.97 

3.68 

3.36 

3.00 

2.61 

2.20 

9 

5.41 

5.32 

5.20 

5.09 

4.96 

4.31 

4.64 

4.45 

4.23 

3.97 

3.68 

3.36 

3.00 

2.61 

2.20 

8 

5.41 

5.31 

5.20 

5.09 

4.96 

4.81 

4.64 

4.45 

4.23 

3.97 

3.68 

3.36 

3.00 

2.61 

2.20 

7 

5.41 

5.31 

5.20 

5.08 

4.96 

4.82 

4.65 

4.45 

4.23 

3.97 

3.63 

3.36 

3.00 

2.61 

2.20 

6 

5.41 

5.31 

5.20 

5.09 

4.96 

4.82 

4.65 

4.45 

4.23 

3.97 

3.68 

3.36 

3.00 

2.61 

2.21 

5 

5.41 

5.31 

5.20 

5.08 

4.96 

4.82 

4.65 

4.43 

4.23 

3.97 

3.68 

3.36 

3.00 

2 . 61 

2.20 

4 

5.41 

5.31 

5.20 

5.09 

4.97 

4.82 

4.66 

4.46 

4.23 

3.98 

3.69 

3.36 

3.00 

2.62 

2.21 

3 

5.41 

5.29 

5.20 

5.09 

4.96 

4.82 

4.66 

4.46 

4.23 

3.98 

3.69 

3.36 

3.00 

2.62 

2.21 

2 

5.41 

5.29 

5.20 

5.09 

4.96 

4.82 

4.66 

4.46 

4.23 

3.98 

3.69 

3.36 

3.00 

2.62 

2.21 

1 

5.41 

5.29 

5.20 

5.09 

4.96 

4.82 

4.66 

4.46 

4.23 

3.98 

3.69 

3.36 

3.00 

2.62 

2.21 

SIGNIFICANT  WAVE  HEIGHTS  IN  METERS 
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20 

1.73 

1.35 

.93 

.52 

.21 

19 

1.78 

1.35 

.93 

.52 

.21 

13 

1.78 

1.35 

.93 

.52 

.21 

17 

1.78 

1.35 

.93 

.52 

.21 

16 

1.78 

1.35 

.92 

.52 

.21 

15 

1.78 

1.35 

.92 

.52 

.21 

14 

1.78 

1.35 

.92 

.52 

.21 

13 

1.78 

1.35 

.92 

.52 

.21 

12 

1.78 

1.35 

.92 

.52 

.21 

11 

1.78 

1.35 

.92 

.52 

.21 

10 

1.78 

1.35 

.92 

.52 

.21 

9 

1.78 

1.35 

.92 

.52 

.21 

8 

1.78 

1.35 

.92 

.52 

.21 

7 

1.78 

1.35 

.92 

.52 

.21 

6 

1.78 

1.35 

.92 

.52 

.21 

5 

1.78 

1.35 

.92 

.52 

.21 

4 

1.78 

1.35 

.93 

.52 

.21 

3 

1.78 

1.35 

.93 

.52 

.21 

2 

1.78 

1.35 

.93 

.52 

.21 

1 

1.78 

1.35 

.93 

.52 

.21 

E5 


Coefficients  used  in  the  formulation  of  the  Phillips  atmospheric  in¬ 
put  term 

Wave  phase  velocity 
Coefficient  of  bottom  friction 
Wave  group  velocity 
Depth  of  permeable  layer 
Spectral  energy  in  { f , © )  space 
Wave  frequency 

Wave  frequency  at  the  spectral  peak 
Spectral  energy  in  k  space 
Acceleration  of  gravity 
Water  depth 
Wave  number  vector 

Ratio  of  wave  frequency  to  wind  speed 
Coordinate  distance  along  a  wave  ray 
Source  or  sink  term 
Time 

Wind  speed 
Water  velocity 
Coordinate  cross-shore 
Coordinate  alongshore 
Spectral  parameter 

Linear  growth  coefficient  in  Phillips'  term  for  atmospheric  input 
Permeability  coefficient 
Spectral  parameter 

Wave  direction  or  difference  between  wave  and  wind  direction 
Constant  =  3. 141 59  - . . 

Density  of  air 
Density  of  water 
Spectral  constant 
Spectral  parameter 
Wave  frequency  =  2*f 


